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Abstract: The Oppenheimer-Snyder model, which describes the gravitational collapse 
of a ball of dust, can be used to model thermalization of strongly coupled systems us¬ 
ing the AdS/CFT-duality. It can be used as an alternative to the previously widely used 
thin-shell model in the context of holographic thermalization. We solve the dynamics 
of the Oppenheimer-Snyder collapse in asymptotically AdS space-times and compute the 
two-point function of a boundary spectator field in the geodesic approximation. As a com¬ 
parison we perform the same computation in the thin-shell model and discover oscillatory 
solutions of the thin-shell model for certain equations of state. 
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1 Introduction 

Holographic methods have been well established in the recent years as a serious tool to 
understand the dynamics and thermalization of strongly coupled systems. The AdS/CFT 
duality [1-3] maps the complex process of thermalization of a strongly coupled Af = 4 SYM 
held theory in four dimensions into the relatively simple process of gravitational collapse 
in 5D gravity. The use of the duality has spawned an abundance of interesting work 
attempting to gain insight to different aspects of thermalization (e.g. see references 8-37 in 
[4]). Work that has attempted to describe states close to those in heavy ion collisions has 
been performed in e.g. [5-9], in the sense that the initial state and the distribution of the 
energy-momentum in the boundary held theory is close to the one in heavy ion collisions. 
Although certain simple quantities, such as the aforementioned energy-momentum tensor, 
are easily computed in these numerical works, the computation of other observables can 
be very difficult due to the numerical nature of the work. For computing other observables 
analytical models of the thermalization process are thus favoured. The general theme of 
these analytical models is then that although their symmetries prohibit them describing the 
initial state of a heavy ion collision accurately, they are simple enough that the gravitational 
dynamics can be understood analytically and computations involving different observables 
can be performed in these backgrounds. 
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The simplest analytical model of gravitational collapse used in holography is the thin- 
shell model [10, 11]. It consists of an infinitely thin shell of matter with some equation of 
state collapsing into a black hole. A variation of this simple model is the Vaidya metric, 
where the shell is null (see [12] and references therein). The relative simplicity of the 
thin-shell model enables the computation of different quantities such as higher correlations 
functions of the stress tensor [13], photon and dilepton production rates [14, 15], the two- 
point function in the geodesic approximation [16, 17], and different kinds of entropies 
[4, 18-20] just to name a few. 

The initial state described by the thin-shell collapse however has a peculiar feature: 
Due to the translational symmetry of the setup, all the one-point functions in the boundary 
held theory have their equilibrium value throughout the evolution. This is in stark contrast 
with the initial state of heavy ion collisions, where the initial state is described by a very 
inhomogeneous distribution of matter. Thus the thin-shell model cannot describe perhaps 
the most important observable, the thermalization of the energy-momentum tensor. Thus 
it is clear that the thin-shell model cannot describe in detail heavy ion collisions or other 
processes encountered in the laboratory setting. Instead, we hope that by studying these 
analytical models, we can learn something qualitative about the dynamics of the strongly 
coupled system. When we compute observables, we hope that their features are generic for 
thermalization processes, and can be applied to other scenarios as well. 

It is however not clear which features of the predictions of the thin-shell model are 
generic and which very specific to that setup. For example, the thin-shell model has an 
energy-momentum tensor that is zero throughout in the bulk except for a delta function 
contribution of the shell, and, in fact, certain observables such as the photon and dilepton 
production rates hves been observed to have oscillatory features [14, 15] for which it is 
unclear whether they are generic signatures of thermalization or whether they are caused 
by the rather singular nature of the shell. For the forementioned reasons it is extremely 
important to employ different analytical models for holographic thermalization. Ideally, 
we could compute observables in different scenarios and thus find out which predictions 
seem generic patterns of thermalization and which features are specific to each model. 

An alternative to the thin-shell model is the Oppenheimer-Snyder (OS) collapse. This 
describes the gravitational collapse of a homogeneous and isotropic ball of dust initially 
at rest. This model was first published already in 1939 [21] as a model of star formation 
in 3+1 dimensional asymptotically Minkowski space-time. In the context of AdS/CFT it 
was first mentioned in the appendix of [22], but aside from the equations of motion of the 
collapse it was not developed further there. The advantage of the Oppenheimer-Snyder 
collapse in the context of holography is that since it describes a ball of dust instead of an 
infinitely thin shell of matter, the energy-momentum tensor does not have a delta function, 
but rather is finite everywhere. Thus the resulting space-time is more regular at the edge 
of the ball than at the location of an infinitely thin shell. The model is still symmetric 
enough that for example the equations of motion of a spectator field in the bulk can solved 
at least for some cases. 

It is the study of the Oppenheimer-Snyder collapse as a model of holographic ther- 
malization that this work is devoted to. We will first introduce the model and solve its 
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gravitational evolution. Then as a test observable, we compute the equal time two-point 
function in the geodesic approximation. To compare the OS setup with the thin-shell 
model, we perform the same computations in the thin-shell model. As a by-product of the 
computation, we discover previously unknown oscillatory solutions of the thin-shell model 
in global AdS. 

The paper is organized as follows: In section 2 we introduce the Oppenheimer-Snyder 
model in asymptotically AdS space-times, and discuss its dynamics. In section 3 we review 
the thin-shell model and review the dynamics of the collapse. In section 4 we compute the 
equal time two-point correlator in the geodesic approximation in both of the models and 
compare the results. Finally in section 5 we conclude and discuss possible developments. 


2 Asymptotically Anti-de Sitter Oppenheimer-Snyder collapse 


The Oppenheimer-Snyder collapse is conceptually the simplest model of gravitational col¬ 
lapse: A homogeneous and isotropic ball of pressureless dust undergoing gravitational col¬ 
lapse. Since we are looking for spherically symmetric solutions, we use global coordinates 
for AdS, where constant-?' hypersurfaces have explicit rotational symmetry. 

Outside the ball the energy-momentum content is zero. The solution with zero X^„, 
negative cosmological constant and rotational symmetry is given by the AdS-Schwarzschild 
solution: 

ds 2 = -f(r)dt 2 + -^—+r 2 dQ 2 n , /(r) = 1 +Ar 2 - , (2.1) 

j[r) r n 1 

where is the metric of a ?r-sphere and A is the cosmological constant. In this coordinate 
system the edge of the sphere has some time dependent value of the radial coordinate r s . 
Thus this coordinate patch is valid for r > r s . 

Inside the ball of dust the simplest coordinate frame to use is the rest frame of the 
dust. Since the dust is homogeneous and isotropic, we want a solution which has only 
time-dependent energy-momentum content, negative cosmological constant and rotational 
symmetry. This metric is of course the hyperbolic Friedmann-Robertson-Walker metric, 


ds 2 


—dr] 2 + a(r/) 2 


' dp 2 

l + P 2 


+ P 2 dnl 


( 2 . 2 ) 


Imposing the Einstein equation with a negative cosmological constant and a pressureless 
fluid at rest and then using the energy-momentum continuity for the fluid, we arrive at the 
Friedmann equation 


da \ 2 2 

dp ) n(n + l)Mpj 


„n+l 

£q 0 1 — A a 2 
a 71 ” 1 


(2.3) 


where eo is the energy density of the dust when the scale factor has the value a q. Since 
this is in the dust’s rest frame, and its total quantity is conserved, the edge of the ball sits 
at a constant value of the radial coordinate p, which we shall call po. Thus the validity of 
this coordinate patch is 0 < p < po- 
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In order to glue these two coordinate patches together, we relate the two different 
patches together at the edge of the sphere and then require that the induced metric there 
is the same in both of the coordinate systems and that the extrinsic curvature is continuous 
across the edge of the sphere. These conditions are the Israel junction conditions [23] with 
zero energy-momentum content at the hypersurface joining the two patches. Although this 
might seem like heavy machinery for the problem at hand, we use it since it is applicable 
later on when we consider the thin-shell case as a comparison to the OS model. 

The coordinates describing the edge of the sphere are chosen to be the angular coor¬ 
dinates, which are identified to be the same inside and outside, and the proper time at the 
edge, t. From the induced metric we can see that we can identify the time coordinate ij at 
the edge of the sphere with the proper time at the edge, i.e. rj s = r. From that it follows 
that r s (r) = poa(r). Additionally the time and radial coordinates t s and r s are related by 

-1 = -f(r s )i 2 s + r 2 /f(r s ) . 


Next we use the Israel junction conditions to relate the curvature on both sides. The 
junction conditions state that the difference of the curvature 1 

K . _ n , r „ 

across the joining hypersurface is related to the energy-momentum content at the hyper¬ 
surface, which in this case is zero. Here X a are the n + 2-dinrensional coordinates of the 
n + 1-dimensional hypersurface parametrized by the n + 1 coordinates called £*, and n a is 
the normal vector of the hypersurface. %j is the induced metric on the surface. Since the 
energy-momentum content is zero, the extrinsic curvatures inside and outisde are related 
by 

(■Kij - K Ifij )inside = ~ ^ 7i? )outside ' (^) 

A straightforward computation gives the non-zero components of the above expression in 
the inside coordinate system to be 


K tt — K 7 tt = — n 


\J 1 + Po 

Poa(r) 


K(j>n<l>n K r Y<t>n(t>n ~ 3<t>n<t>n{ n 1 ) 


n/ 1 + Po 


P 0 ci(t) 

while in the outside coordinate system the non-zero components are given by 


( 2 . 6 ) 

(2.7) 


K tt — K 7 tt = — n 


Vf( r s) + rl 


K 1<Pn<t>n ~ 9<j>n<t>r, 


1 d 

r. r, dr 


Wf(r S ) + r 2 s 


( 2 . 8 ) 

(2.9) 


1 Note that there is an ambiguity in the literature related to the sign of Kij. For example [24] defines 

Kij with an opposite sign. In this case the sign ambguity does not matter, but can cause confusion when 
the energy-momentum of the joining hypersurface is non-zero. 
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Figure 1. A schematic of the coordinate system (a, p) inside the sphere. The grey dashed lines are 
null geodesics. The grey region is the lower half of the manifold, where the sphere is expanding, 
however, the a-coordinate covers only the upper (or lower, depending on the sign choice) region. 
Note that there are null rays emanating from the center of the sphere which hit the singularity 
before encountering the edge of the sphere. 


Equating the expressions of the inside and outside coordinate system, we get the simple 
identity 

f(r s ) + r 2 s = 1 + pi . (2.10) 

Since po is a constant, this is an equation of motion for the radial position of the shell in 
the outside coordinate system, r s . Substituting r s = poa{r) into this equation, we get an 
equation of motion for the scale factor a. If one substitutes the explicit functional form of 
/, this can be seen to be exactly the Friedmann equation, already given in equation (2.3). 

In principle one can proceed with the above derivation with also non-zero pressure of 
the fluid the sphere is composed of, however, one will then see that a solution satisfying 
both equations (2.3) and (2.10) can only be found if p = 0. This is not surprising, since if 
the pressure of the fluid would be non-zero, there would be a pressure gradient at the edge 
of the sphere and the fluid would no longer stay in rest in the inside coordinate system. 

2.1 The motion of the Oppenheimer-Snyder sphere 

Next we analyse the motion and the causal structure of the Oppenheimer-Snyder collapse. 
To do that, we use inside the scale factor as the time coordinate, that is, we make the 


- 5 - 







( 2 . 11 ) 


coordinate substitution da = ^ dr], so that the metric becomes 

2 _ da 2 a 2 2 . 2 2 ;n 2 

as — — + ——+ a P dil n , 

h(a) 1 + p 

where the metric function h(a) = comes directly from the Friedmann equation (2.3). 

It is convenient to relate this function to / via 

h{a) = 1 + 1 ~ ■ (2.12) 

Po 


The time evolution of the diameter of the sphere can be found from equation (2.10). 
The first observation is that since f(r) grows without bound, for any po there is always a 
turning point r max , where f s = 0. Furthermore, since 


d 2 r s 1 d , 2 
dr 2 2 dr s s 


-firs), 


(2.13) 


the acceleration at the turning point is always negative, and at the turning point the sphere 
starts to collapse. This means that there is also a maximum value of the scale factor, a max , 
which is the value of the scale factor at the turning point. The space-time also has a 
space-like singularity at a = 0, which is the “big crunch” of the space-time. 

Due to the time reversal symmetry of the setup, in principle the evolution of the sphere 
is such that it emerges from a singularity and expands until it reaches r max , after which it 
collapses back to a singularity. Thus the coordinate system (a, p) covers only one half of the 
manifold. We consider only the upper half of the space-time which describes the collapse of 
the shell starting form the turning point, with the view that the system is initially prepared 
to the state described by the sphere at the turning point. 

It is also important to note that while there is some value of the “time coordinate” 
a = ah when the horizon is formed, there is no qualitative change inside. In fact, although 
the value of and thus a/,, can be changed by tuning the value of po, the Friedmann 
equation is indifferent to the choice of that parameter. 


3 The thin-shell model 

The thin-shell model has been discussed widely in the literature [4, 13-15, 18, 19] in the 
context of holographic thermalization. Another variation on the same theme is the Vaidya 
metric, where the collapsing shell is null (see e.g. [12]). However, in most of these works the 
collapse has been studied in the Poincare patch of the AdS space. In this coordinate system 
the shell, which is located at constant value of the r-coordinate, is actually an infinite flat 
sheet with translational symmetry. To enable a comparison with the Oppenheimer-Snyder 
collapse, we need to use the global coordinate system of the AdS-space, where a constant- 
r-hypersurface is an actual spherical shell with rotational symmetry. It turns out that the 
spherically symmetric case has interesting behaviour not present in the previously studied 
case in the Poincare patch. The notation follows closely that of [4], so we suggest the reader 
to refer to that one for more details. 
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The metric of the thin-shell model in global AdS is given by 


ds 2 


-f(r)dt 2 + 


dr 2 

f(r) 


+ r 2 dn 2 n 


where the metric function / is now defined by 


(3.1) 


/ 


/_ = 1 + Ar 2 , r < r s 

f+ = 1 + Ar 2 - ^ , r>r s 


(3.2) 


and the location of the shell r s is time dependent. The radial coordinate r and the angular 
coordinates are continuous at the shell, but the time coordinate is discontinuous, and thus 
there are two different time coordinates, t— and t+, the time inside and outside of the shell 
respectively. To parametrize the shell, we use the angular coordinates and the proper time 
of the shell, r, and use an overdot to denote a derivative with respect to r. 

Unlike in the OS collapse, now both the inside and outside solutions are vacuum 
solutions, and all the energy-momentum content is located at the shell. This means that 
the Israel junction conditions are now 


{Kij Klij )inside (K-ij ^ Tij) ou tside — 8ng 5 Sij , (3.3) 

where S l3 is now the energy-momentum tensor of the shell in the ^-coordinate system. We 
use the ideal fluid form, 


S'*- 7 = (p + p)u l v? + py*- 7 , p = cp , 


(3.4) 


where u l = Sq, since the fluid is at rest in the ^-coordinate system. The expressions from 
equations (2.8) and (2.9) now apply both inside and outside, /_ and /+ substituted for 
/ respectively. Plugging these into the junction conditions of equation (3.3), we get two 
independent equations: 


1 d 


9 (fin (fin 


r s r s dr 


“ (v7-(r s ) + r 2 s - s/f+(r a ) + r?) 
rr 1 (y/f-(r s ) +r 2 s - y/f+(r s ) + f|) 


-8ng 5 p 

—8ixg<og§ n §nP 


(3.5) 

(3.6) 


These equations hold for arbitrary relation between p and p. We however assume a constant 
equation of state c, and thus the above equations are solved by 


Vf- + r 2 s - Vf+ + r 2 s 


Mr l—n(i+c) ; 


(3.7) 


where M is now an integration constant characterizing the energy of the shell. From this 
we can also solve the equation of motion for the shell, 


^2 = ^!_2(l-n(l+c)) 
1 s 4 ' s 


/- + /+ r, 2(n(1+c) - 1} (/--/ + ) 2 

2 + 4M 2 


(3.8) 


The integration constant M encodes information about the initial conditions of the shell’s 
motion. Since we want to compare with the motion in the Oppenheimer-Snyder collapse, 
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we assume the shell starts to fall from a turning point ro at r = 0, where f s = 0. We can 
evaluate equation (3.7) at the turning point, and thus we get 


M = 


Vf-{r o) - yY+M 


r 


1—n(l+c) 
0 


This expression with equation (3.8) determines the motion of the thin shell. 


(3.9) 


3.1 Trajectories of the thin shell 

Next we wish to inspect the possible trajectories of the shell. Instead of solving the position 
of the shell r s explicitly as a function of the coordinate time t s , it is sufficient to look at 
the expression r 2 in equation (3.8). First of all, since f 2 s is a square of the physical velocity, 
it has to be positive. Additionally, since r s = \^r'r 2 , the same expression also encodes 
information of the acceleration for the trajectory. Thus by computing the acceleration at 
the turning point, 


d 2 r. 


dr 2 


r s =r o 


mn — l n(l + c) — 1 M 2 
r oA + -w A I - 


n c m r. 


2 J2cn—1 


0 


n(l+c) 

r o 


4 M 2 


(3.10) 


we can determine towards which direction the shell starts to move. As long as the expression 
in the square brackets is positive, the shell starts to fall at the turning point. However, if 


n c mn r, 


2 Jlcn— 1 


0 


M 2 


> 4roA + (n - 1)^ + [n( 1 + c) - 1] 


'o 


n(l+c) 

r o 


(3.11) 


instead of collapsing, the shell starts to expand. Thus all trajectories can be divided to 
two categories: The ones where the shell starts to collapse, and the ones where the shell 
starts to expand at its turning point ro- 

Next question then is, what happens to the shell after its initial turning point? Is there 
another turning point, or will the shell always either collapse to a black hole or expand all 
the way to the boundary? The EOM cannot be solved analytically for turning points for 
general values of n and c. When investigating the solutions numerically, it turns out that 
both kinds of behaviour are present. Thus the solutions can be divided into two additional 
categories: The ones which have only a single turning point, and the ones which have more. 
Thus we have four general categories of solutions, illustrated in figure 2: 


1. Shell collapses to a black hole (solid line) 

2. Shell starts to collapse, but has a turning point at smaller radius (dotted line) 

3. Shell starts to expand, but has a turning point at a greater radius (small dashed line) 

4. Shell expands all the way to the boundary (large dashed line) 

In case 1, the expression in equation (3.8) has only a single root. For case 2, the expression 
has three roots, however, since the initial condition is given at ro, the smallest root is never 
reached, and the shell ends up oscillating between the ro and the middle root. For case 3 
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Figure 2. The speed squared of the thin shell for different equations of state. Solid line c = 0.328, 
the dotted line c = 0.33, the small dashed line c = 0.3323 and the large dashed line c = 0.3334. 
The divergence of the velocity at r h is a sign of a coordinate singularity at the horizon. The other 
parameters have been chosen to be m = 5, A = 1, ro = 10. 

the expression has again three roots, and now the shell oscillates between r q and the largest 
root. For case 4, the the expression has two roots, but the smaller one is never reached, 
and instead the shell starts and continues to expand from ro towards infinity. 

The different trajectories in figure 2 differ only by equation of state. In the case of 
c > — 1 and n = 3, one can look at the powers of the different terms in equation (3.8): At 
small r s the first term will dominate no matter what the value of c. However, at large r s , 
if c < ^ the middle term will dominate and the expression will turn negative at sufficiently 
large r s , and if c > ^ then the last term will dominate and the expression will be positive 
at large r s . Thus one critical value of c is which separates cases 3 and 4. The critical 
values of c separating cases 1 and 2, and 2 and 3, respectively, are however dependent on 
the other choices of parameters and initial conditions, and cannot be solved analytically. 
Thus the exact numerical values of the choices of c in figure 2 are not important, but 
instead it is sufficient to note that if one denotes c* as the equation of state of the different 
cases, then ci < C 2 < C 3 < C 4 , C 3 < ^ and C 4 > K 

The trajectories of the shell encode three qualitatvely different effects: Geodesic motion 
in AdS, gravitational pull of the shell, and the pressure of the shell. The non-linear interplay 
of these different effects is the cause of this rich structure for the possible trajectories. The 
collapsing shell (case 1) and the eternally expanding shell (case 4) are the most familiar 
solutions, where as the eternally oscillating shell (cases 2 and 3) are more peculiar. In 
the context of holography it would be tempting to try to interpret these solutions as non- 
thermalizing states in the boundary field theory. However, when c / 0, the shell is made 
out of matter that has pressure only in some spatial directions, and not in others. The 
classical analogy of this would be a sheet of rubber or a network of springs. In order for 
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this shell to have a holographic interpretation in the boundary field theory, there must be 
a bulk degree of freedom from which the shell is constructed. Whether this is possible 
for the non-zero values of c is an unanswered question. One familiar case where there is 
pressure along the shell is an instanton-like solution of e.g. a scalar field in the bulk, but 
for these solutions the equation of state parameter is usually c ~ — 1. 

As a final note, we must note that for the equation of state c = 0, which one would 
imagine to be easiest to construct from the degrees of freedom available in the holographic 
description, the shell always collapses. Due to the lack of pressure, this is also the configu¬ 
ration we will use when making comparisons with the Oppenheimer-Snyder collapse in the 
following sections. 

4 The equal time two-point correlator 

To probe the thermalization described by the Oppenheimer-Snyder collapse and the thin- 
shell collapse, we compute the equal time two-point correlator of an operator with a large 
conformal weight using the geodesic approximation. The geodesic approximation (for de¬ 
tails and derivation see e.g. [25, 26]) is essentially the saddle point approximation in the 
limit of large scaling dimension of the operator whose correlator we are computing. The 
description for the correlator is then given by 

(O(x)O(x')) ~ e- 2A e- Af7 ( x ’ x ') (4.1) 

where A is the conformal weight of the bulk particle or equivalently the scaling dimension 
of the dual operator, cr(x,x / ) is the length of the geodesic connecting the two boundary 
points x and x 7 , and e is a cutoff parameter regularizing the length of the geodesic as it 
approaches the boundary. From this we see that the computation of the equal time two- 
point correlator is reduced to computing the geodesic length between two points residing 
in the boundary at equal times. Due to the spherical symmetry of the setup, the geodesic 
length, and thus the correlator, will only depend on the time coordinate of the points and 
one angle separating the points in the n-sphere which forms the boundary. This angle we 
shall call <p. 

The cutoff e corresponds to some (large) finite radius R to which the geodesics are 
computed. Since the divergence of the geodesic length when taking R —> oo is present even 
in the case of empty AdS, we define the renormalized geodesic length by 

2 

a(t, ip) = a(t, p) -— In2\ // A R, (4.2) 

VA 

where the last term on the RHS is the divergent part of the length of a spacelike geodesic 
in empty AdS. The renormalized a is thus finite in the limit R —> oo. The renormalization 
scheme here is arbitrary, and in fact a will depend on the choice of the scheme. However, 
we will only compare the lengths of geodesics with matching boundary endpoints, and the 
choice of the renormalization scheme does not affect the difference of the lengths of these 
geodesics. Similar scheme was used in e.g. [12], 
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The derivation of the matching conditions of the geodesics in this section follows the 
same recipe and computational machinery described in [4] , appendices of which the reader 
is referred to for more detailed derivations. 


4.1 Calculational strategy 

We want to know the length of the shortest space-like geodesic connecting two boundary 
points with the time coordinate t and the angular separation ip. Although a simple state¬ 
ment, this quantity is not trivial to compute, since if we start from the boundary points we 
do not know to which direction the geodesic should start to propagate from the boundary 
- that behaviour encodes the bulk geometry of the space-time. Instead, we approach the 
problem from a different point of view: A geodesic that connects the two points with the 
same time coordinate must have a turning point somewhere in the bulk, where both the 
radial and temporal velocities 2 are zero. Thus we can start from any point in the bulk 
and look for geodesics which have a turning point at that point, and then propagate these 
geodesics all the way boundary. This means that for each set of initial conditions (in the 
Oppenheimer-Snyder case (a,p) and in the thin shell case (i_,f)) we can compute the 
time t when the geodesic hits the boundary, the angular separation ip, and the length of 
the geodesic u. If we then scan through all possible values of the initial turning point, we 
should recover a geodesic length for any possible combination of (p and t. 

It turns out that one needs to consider also geodesics which penetrate the edge of 
the sphere or the thin shell after the horizon has formed, i.e. geodesics which cross the 
coordinate singularity at r = in the outside coordinate patch. Although the coordinate 
system has a singularity at that point, with careful regularization all physical quantities 
will stay finite and the computation can be carried out in this coordinate system. A more 
robust approach is to use the Eddington-Finkelstein coordinate system outside, details of 
which are give in Appendix B. 


4.2 Geodesics in the Oppenheimer-Snyder collapse 

To calculate the length of geodesics connecting two equal time boundary points in the 
Oppenheimer-Snyder geometry, we need to solve the set of geodesic equations in the two 
coordinate systems defined by equations ( 2 . 1 ) and ( 2 . 11 ), and then match the solutions at 
the edge of the sphere of dust. 

The geodesic equations in the inside coordinate system are given by 


h! . o o - 9 ah p 2 
-—a +ap hip + —— 9 + a 
2 h 1 + p- 


-p(i + p 2 )4> 2 + 


2 a p p p 2 
a 1 + p 2 
2 a 4> 2 pf> 


P = 


+ + <f> 

a p 


0 

0 

0 , 


(4.3) 

(4.4) 

(4.5) 


2 Here we use the term velocity to denote a derivative of the coordinates of the geodesic with respect to 
the proper length of the geodesic. Although these quantities are not velocities in the physical sense, they 
appear in the geodesic equations in the same fashion as actual velocities would for time-like geodesics. 
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Figure 3. Comparison between the trajectories of the pressureless (c = 0) thin shell and the OS 
collapse with equal initial radius ro as a function of the proper time of the shell. Solid line is the 
position of the edge of the dust sphere of the OS collapse while the dashed line is the location 
of the thin shell. The dotted line is the location of the horizon. The difference between the two 
trajectories gets even smaller if rh is decreased or ro increased. Here the other parameters A = 1, 
m = 10, n = 3 and rg = 2 have been chosen such that the difference is big. 


where ' denotes a derivative w.r.t. a and " a derivative w.r.t. the affine parameter of the 
geodesic, which we choose to be the geodesic length a. 

The geodesic equations (4. 3-4. 5) can be integrated to give 


L 


a 2 p 2 


a = 


P = 


h(a) 

1 + p 2 r 


E 

-2 - 1 
a z 


E 

P 


(4.6) 

(4.7) 

(4.8) 


where E and L are integration constants, i.e. constants of motion. Their names are chosen 
to remind the reader of energy and angular momentum, since they are the analogues of these 
quantities. The values of the constants of motion are determined from the requirement that 
the geodesic has a turning point at (a,p): 


a 2 (a = a) = 0 =>- E = a 2 

p 2 (p = p) = 0 => L = ±ap 


(4.9) 

(4.10) 


Now that we know how the geodesic propagates from the turning point to the edge of 
the sphere, we need to solve how the geodesic propagates outside the sphere, and how the 
geodesic solution inside is matched to a solution outside. The geodesic equations outside 
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are given by 


f'rt 

f 


+ t = 0 


f -2 , 1 f ft 12 


~Yf*" Z + f t z ~ r f 4> 2 + r = 0 


2 r<j) 


+ <t> = 0 . 


These can be integrated to give 


fi = E 
r 2 <j) = L 



(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 


To relate the constants of motion outside to the trajectory of the geodesic inside, we 
construct a coordinate system (r, A) which is continuous at the edge of the sphere and 
its close surroundings and use that to transform the geodesic derivatives in the inside 
coordinate system to those in the outside coordinate system. This computation is described 
in more detail in Appendix A, but its result is given by 


E = 





L = L = ±ap . 


(4.17) 

(4.18) 


The procedure to compute the geodesic length for all times and angles is then as follows: 
Given the initial conditions a and p, integrate equations (4.6 - 4.8) until the geodesic 
reaches the edge of the sphere. Then use the matching conditions (4.17) and (4.18) to get 
the constants of motion outside. With the constants of motion, integrate equations (4.14 
- 4.16) until the geodesic reaches the boundary, which is represented by some (large but 
finite) cutoff in the radial coordinate, R. Repeat the computation for all possible values of 
the initial conditions a and p. One must discard those cases where either the geodesic does 
not exit the sphere before a = 0, or the geodesic does exit the sphere, but turns around 
and falls back to the sphere and thus never reaches the boundary. 


4.3 Computation in the thin-shell model 

The computation in the thin-shell case is very similar to the Oppenheimer-Snyder case. 
Since the metric in equations (3.1) and (3.2) is the same as the metric of equation (2.1), 
the geodesic equations both inside and outside are given by the equations (4.14 - 4.16). We 
shoot the geodesic from its turning point (t_,r_). From the requirement that r = 0 and 
t = 0 at the turning point, we are able to derive that the constants of motion inside are 
given by 

= 0 and L_ = ±f . (4-19) 
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We can then integrate the geodesic equations (4.14 - 4.16) until the shell, r = r s , at which 
point we need to switch to the geodesic equations in the outside coordinate patch. These 
are of course given by equations (4.14 - 4.16) as well, but now with f = f + and with 
different constants of motion. These constants can be derived (for more detail, see e.g. the 
Appendix B of [4]) to be 


E + = 




L+ = r . 


(4.20) 

(4.21) 


The procedure to compute the geodesic length is after this point very similar to the 
Oppenheimer-Snyder case: Given the initial radial coordinate f, integrate the equations 
for some value of r s . This corresponds to a choice of initial time. Then use the matching 
conditions (4.20-4.21) to compute the constants of motion outside, after which integrate 
the geodesic equations (4.14-4.16) until the geodesic reaches the ’boundary’, i.e. the ra¬ 
dial cutoff R. Then rinse and repeat for all values of f and r s . Again, one must discard 
geodesics which never reach the boundary and also take into account geodesics that exit 
the shell after the horizon has formed. 


4.4 Results 

We wish to compare the thermalization of Oppenheimer-Snyder to that of a thin-shell 
collapse. To make the comparison of the two qualitatively different models meaningful, we 
make the mass of the final black hole and the initial radius of the sphere and shell to be 
the same for both models. This corresponds to setting the final temperature of the state of 
the boundary field theory as well as the angular scale of the initial state to be the same for 
both states. Also, since the sphere of dust is pressureless by construction, we also choose 
the shell to be pressureless, i.e. c = 0. The differences in the correlator for the two different 
models are then due to two distinct sources: different speed of the gravitational collapse 
and different trajectories of geodesics due to the different geometries. 

The motion of the radius of the dust sphere and the shell are compared in figure 3. 
Both are released at t = 0 and start to collapse. The Oppenheimer-Snyder collapse is 
always faster than the pressureless thin-shell collapse. This can be understood by the 
distribution of the gravitating matter, since in the OS there is matter in the middle which 
attracts the outer layers as well. The location of the shell or the edge of the sphere 
corresponds to the smallest equilibriated scale in the boundary theory, and thus a given 
scale always thermalizes first in the Oppenheimer-Snyder collapse. It should be noted that 
the parameters in figure 3 have been chosen to make the difference in the evolution of the 
two different models to be large, and for generic choices of parameters the actual difference 
between the two trajectories is even smaller. 

Even if the difference between the two different trajectories is small, it does not im¬ 
ply that the difference in the geodesic length between two boundary points is small as 
well. In the thin-shell model, all the energy-momentum resides at one radial location, r s , 
which causes a jump in extrinsic curvature normal to the shell. This jump in curvature 
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Figure 4. Renormalized geodesic length as a function of the angular separation of the boundary 
points evaluated at different times. Solid line is the equilibrium value, dotted the thin-shell collapse 
and dashed the Oppenheimer-Snyder collapse. From bottom to upwards the lengths are evaluated 
at times t = 0, 0.5 and 0.75. The values of other parameters are n = 3, m = 10, rp = 3 and A = 1. 



0 0.2 0.4 0.6 0.8 

t 


Figure 5. The unnormalized correlator e~ G for boundary angular separation <f> = 1.76 as a function 
of time. The solid line corresponds to the equilibrium case, the dashed to the Oppenheimer-Snyder 
case and the dotted to the thin-shell collapse. Other parameters are the same as in figure 4. 


causes the geodesic to refract as it crosses the shell. In the Oppenheimer-Snyder model 
there is no delta function in the energy-momentum tensor and correspondingly no jump in 
extrinsic curvature which would cause the geodesic to refract. Thus the geodesics in the 
Oppenheimer-Snyder model are more smooth, and geodesics that penetrate the shell devi¬ 
ate less from the equilibrium then in the thin-shell case. It should be noted that due to the 
coordinate systems in use, this smoothness is not apparent when numerically computing 
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Figure 6. The unnormalized correlator e -cr for boundary angular separation cf> = 2.94 as a function 
of time. The solid line corresponds to the equilibrium case, the dashed to the Oppenheinrer-Snyder 
case and the dotted to the thin-shell collapse. Other parameters are the same as in figure 4. 



the trajectories of the geodesics. 

In figure 4 we have plotted the geodesic length between two points as a function of 
the angular separation of these two points. The solid line is the equilibrium value, the 
dashed lines the Oppenheimer-Snyder model and the dotted lines the thin-shell mode. The 
pairs of dashed and dotted lines correspond to different boundary times. First one can 
observe that the point where the dashed and dotted line meet the solid line corresponds 
to the largest scale that has equilibriated, which corresponds to the shell motion in figure 
3. Furthermore, the Oppenheimer-Snyder geodesics are always closer to equilibrium for all 
angular scales compared to the thin-shell case. Finally, one can note that the difference 
between the two models is largest in intermediate times, and very early and very late in 
the gravitational collapse the difference is smaller. 

In figures 5 and 6 we have plotted the exponential of the geodesic lengths of a fixed 
angular separation as a function of time. According to equation 4.1, this is directly pro¬ 
portional to the equal time two-point correlator, modulo an overall normalization factor 
that has not been determined. Again one can see that a given scale thermalizes faster in 
the Oppenheimer-Snyder model than in the thin-shell case. It should be noted that the 
thermalization starts with zero velocity, corresponding to the shell or sphere being initially 
at rest, but reaches the thermal value in finite time of the boundary observer. 

When computing the geodesic length, for intermediate times and large angular sepa¬ 
ration (that is, close to p = n) we found the geodesic length to be a multivalued function. 
That is, it appears that for a given angular separation cp, there are several geodesics that 
connect these two boundary points. To illustrate this we have plotted the lengths of all 
geodesics connecting boundary points with the angular separation cf = 2.88 in figure 7. One 
can see that for intermediate times there are three different lengths corresponding to this 
fixed angular separation. If the geodesic approximation is taken at face value, the shortest 
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Figure 7. The geodesic length between boundary points with angular separation </> = 2.88 as a 
function of time in the OS collapse. For intermediate times there are several geodesics with different 
lengths connecting the boundary points. The parameters are chosen to be A = 1, m = 4, n = 3 
and ro = 4. 


geodesic is the dominant one, and thus the two-point correlator becomes discontinuous. 
This feature is present in both the Oppenheimer-Snyder and the thin-shell models. The 
region where the geodesic length seems to be multivalued is always limited to intermediate 
times and large angles, i.e. angles close to n. Similar behaviour was found in the Vaidya 
model in [12], Specifially figure 7 should be compared with figure 7 of [12]. There the 
authors argue that the geodesic approximation breaks down, and the multivaluedness of 
a is not to be taken seriously. In the present work we have followed the authors of [12] 
and assumed that the value of the correlator is given by the geodesic length given by the 
largest values in figure 7 and thus it is continuous. Indeed, when we have produced the 
plots, we have explicitly selected the points corresponding to the continuous values, and 
discarded the additional points. 

Although we can make the correlator thus continuous, its derivative is not: By looking 
at figure 7 it is clear that if we choose the two curves that extend to small and large times, 
their intersection (corresponding to f ~ 5.5 in figure 7) has a kink. Thus for small angles 
where the geodesic lenght is not multivalued, the correlator is smooth as evident in figure 
5, but for large angles, it is not, as demonstrated by the kink in figure 6. The kink is 
present in both the Oppenheimer-Snyder case (dashed line) and the thin-shell case (dotted 
line). 

The discontinuity of the derivative of the correlator has been discussed already in 
e.g. [17], where it was argued to be physical. It is interesting to note, however, that both 
in [17] and in [12] the computation was carried out in the Vaidya metric, which has a 
null hypersurface separating the inside and outside space-times, i.e. the matter collapses 
with the speed of light. In this work, however, we have not only looked at time-like shells 
collapsing, but also the Oppenheimer-Snyder collapse where there is no shell, and thus 
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no discontinuity in the curvature and no delta function inthe energy-momentum tensor. 
Nevertheless the jump in the time derivative of the correlator is present, demonstrating 
that the jump does not stem from the null surface in Vaidya nor from the delta-function 
like distribution of matter. If one assumes that the physical correlator should be smooth, 
then it seems that the geodesic approximation breaks down not only in the sense that it 
produces multivaluedness, but also since it does not produce a smooth correlator. 

5 Conclusions 

The Oppenheimer-Snyder model describes the gravitational collapse of a ball of dust and 
can be used as an analytical model for holographic thermalization. It serves as an alterna¬ 
tive to the previously commonly used thin-shell models. One advantage of the OS model is 
that the OS metric is smoother than the thin-shell model metric due to the lack of a delta- 
function in the energy-momentum tensor. More important however is that the OS model 
offers a different analytical model for holographic thermalization: As we are interested in 
the general features of thermalization in holographic models and not features specific to 
the thin-shell model, the OS model can function as a new comparison to the thin-shell 
model, and we can hope that the common features of all these models are generic and can 
be applied to the physical applications of strongly coupled systems. 

We investigated the speed of the collapse in the OS model compared to the thin-shell 
case. We found that the OS model collapses always faster than the corresponding thin-shell 
model due to the more homogeneous distribution of the energy-momentum. The speed of 
the collapse is however not an observable, and thus we also computed the equal time two- 
point function in the geodesic approximation for both cases. With similar initial conditions 
(same final mass, same initial radius), the OS correlator was found to be at all times closer 
to equilibrium than the thin-shell correlator. This was due to two separate aspects: First of 
all the OS model collapses more quickly, but furthermore, since the metric is more smooth 
at the edge of the sphere than at the location of an infinitely thin-shell, the geodesics are 
more smooth and refract less, and thus deviate less from the equilibrium geodesics. Thus 
we concluded that the OS model always thermalizes faster then the thin-shell model. 

When we computed the evolution of the radius of the thin shell to compare with the 
OS model, we found that for certain non-zero equations of state the thin shell in global 
AdS has curious trajectories: In addition to the familiar behaviour of a collapse to a black 
hole or an infinite expansion all the way to the boundary, we found that there are solutions 
where the shell ends up oscillating between two turning points, thus never collapsing to a 
black hole but never dissipating to inihnity either. It is unclear what the interpretation of 
these states are in the boundary held theory, or if these cases correspond to such equations 
of state which are unattainable in realistic applications of the duality. 

When computing the two-point function, we also found that for large times and large 
angular separations, the geodesic length is multivalued and/or discontinuous. This feature 
of the geodesic approximation has been previously found in the literature in the Vaidya 
model [12], but here we demonstrated for the first time that this discontinuity exists even 
in the metric which has no delta-functions in the energy-momentum tensor. It seems that 
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the geodesic approximation of the two-point correlator breaks down: not only does the 
correlator seem multi-valued and discontinuous, but there is also the qualitative issue that 
the correlator seems to know about physics inside the horizon and thus violating naive 
understanding of causality. Thus it would be interesting to investigate where exactly the 
approximation breaks down. If we were investigating the retarded correlator, one could 
solve a wave equation with certain boundary conditions to obtain the full correlator and 
then see how this compares with the approximation. However, for the correlator at hand 
such a simple calculation is not possible. Instead, a procedure described in e.g. [27] is 
required, which given the non-trivial metric seems a very difficult task. 

The main interest in this and previous work is to discover which features of the an¬ 
alytical models are generic enough that they can be expected to be present in physically 
realistic scenarios (e.g. heavy ion collisions). For the equal time two-point correlator we 
found out that results are qualitatively very similar in both OS and thin-shell models. We 
also found that the apparent breakdown of the geodesic approximation is present in both 
models. The next step would be to compute the different observables for the OS model 
which have already been computed in the thin-shell model in the literature [13-15] and see 
which features of the thermalization are in fact just features of the thin-shell model. For 
example the photon production rate and dilepton spectrum computed in [14, 15] exhibited 
oscillations which were speculated to be related to the presence of the delta-function-like 
shell. Performing the same computation in the OS model would immediately determine 
the source of the oscillatory behaviour. 
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A The continuous coordinate system 

To compute the matching conditions, e.g. how a geodesic travers the shell from one coor¬ 
dinate system to another, we construct a coordinate system which is explicitly continuous 
at the edge of the sphere. This coordinate system uses coordinates r and A as its time 
and radial coordinates. A is defined as the proper distance from the hypersurface which 
is the edge of the sphere normal to the hypersurface. In more practical terms, it is de¬ 
fined as the length of the geodesic that at the hypersurface is normal to the hypersurface, 
i.e. dx^/dX = n The time coordinate r is then defined to be the proper time of the shell 
when that that geodesic intercepts the hypersurface. Although this coordinate system has 
a very complicated form in the bulk, at the hypersurface joining the inside and outside 
metric it takes the form of a Minkowski metric. We can then transform the geodesic both 
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inside and outside to this coordinate system and propagate the geodesic infinitesimally 
across the hypersurface. Since the metric is continuous, the geodesic evolves smoothly. 
The only ingredients we then need for this computation are the partial derivatives of the 
different coordinates, which we can use to perform this transformation. For the outside 
(and inside in the thin shell case) metric these are given by 


For the FRW metric inside the sphere in the Oppenheimer-Snyder case, the derivatives are 
given by 


Z =0, ^ =0 and ^ = 


B Eddington-Finkelstein coordinate system 

When computing the geodesics for the equal time two-point correlators, one needs to solve 
for geodesics that cross the shell or leave the sphere only after the horizon has been formed. 
Easiest way to avoid the divergences associated with the singularity at the horizon is to 
use the Eddington-Finkelstein metric for the outside metric, given by 

ds 2 = —f(r) dv 2 + 2 dv dr + r 2 dQ 2 . (B.l) 

The geodesic equations for a space-like geodesic in this coordinate system can be integrated 
to be 

r 2 d) = L (B.2) 


r 2 = f{r) \l-^+E 
E + r 


The location of the shell (v s (r), r s (r ),...) obeys 


fs =n(fVs-T 


and v q = 


r s + V / + r 2 s 


The normal vector of the shell is 


n r = v s f - — 
V v x 


V 

n = Vo 


The coordinate transformation to the continuous coordinate system (r, A) at the shell or 
the edge of the sphere is then given by 




r s + \/f + r 2 s 
f 

r s + Vf + r 2 s 
f 
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